Are braneworlds born isotropic? 
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It has recently been suggested that an isotropic singularity may be a generic feature of brane 
cosmologies, even in the inhomogeneous case. Using the covariant and gauge - invariant approach we 
present a detailed analysis of linear perturbations of the isotropic model Tb which is a past attractor 
in the phase space of homogeneous Bianchi models on the brane. We find that for matter with an 
equation of state parameter 7 > 1, the dimensionless variables representing generic anisotropic 
and inhomogeneous perturbations decay as t — > 0, showing that the model Tb is asymptotically 
stable in the past. We conclude that brane universes are born with isotropy naturally built-in, 
contrary to standard cosmology. The observed large - scale homogeneity and isotropy of the universe 
can therefore be explained as a consequence of the initial conditions if the brane -world paradigm 
represents a description of the very early universe. 

PACS numbers: 98.80.Cq 



INTRODUCTION 

A classic problem of cosmology is finding ways to ex- 
plain the very high degree of isotropy observed in the Cos- 
mic Microwave Background (CMB). In general relativity, 
where isotropy is a special rather than generic feature of 
cosmological models, we need a dynamical mechanism 
which is able to produce isotropy. One of the most effi- 
cient mechanisms for isotropising the universe is Inflation 
[3fj| but it requires sufficiently homogeneous initial data 
in order to start at all 0- Although one can adopt the 
view that it is sufficient to have one such homogeneous 
enough patch in an otherwise inhomogeneous initial uni- 
verse to explain what we observe, this seems somehow 
unsatisfactory. In our view the isotropy problem remains 
open to debate in standard cosmology. 

Over the past few years the brane- world paradigm has 
received considerable attention as apossible candidate 
for string inspired cosmology (see [3] for a recent re- 
view). In this scenario the observable universe is a 4- 
dimensional slice, the brane, embedded in a higher di- 
mensional spacetime called the bulk. Here we consider 
the formulation developed in Q in order to generalise a 
previous model by Randall and Sundrum where the 
bulk is 5 - dimensional and contains only a cosmological 
constant, assumed to be negative. 

In a series of rece nt pape rs a number of authors 
llSSSHmHlHElOil have presented a de- 
tailed description of the dynamics of homogeneous and 
anisotropic brane- worlds, finding a remarkable result: 
unlike standard general relativity, where the generic cos- 
mological singularity is anisotropic, the past attractor 
for homogeneous anisotropic models in the brane is a 
simple Robertson - Walker model T\,. More significantly, 
in I9I llfl it was shown that this result holds true in 



Bianchi IX models, as well as for some simple inhomo- 
geneous models. Since the Belinski - Lifshitz - Kalatnikov 
(BLK) conjecture ^(| suggests that the Bianchi IX be- 
haviour in the vicinity of the singularity is general, i.e. 
that the approach to the cosmological singularity in a 
generic inhomogeneous universe model should locally be 
the same as in Bianchi IX, it has been suggested that the 
isotropic singularity could be a generic feature of brane 
cosmological models. 

If this conjecture 0,0] could be proved correct, brane 
cosmology would have the very attractive feature of hav- 
ing isotropy built in. Inflation in this context would still 
be the most efficient way of producing the fluctuations 
seen in the CMB, but there would be no need of special 
initial conditions for it to start |l7| . Also the Penrose 
conjecture 18] on gravitational entropy and an initially 
vanishin g W eyl tensor would be satisfied in these models 
(c.f. [11123). 

In this paper we prove that this conjecture is true, 
within a perturbative approach and in the large-scale 
and high energy regime, as justified below. We arrive 
at this result through a detailed analysis of generic lin- 
ear inhomogeneous and anisotropic perturbations I2HE2I 
Eii E4I ] of this past attractor T b - This is done by us- 
ing the full set of linear 1+3 covariant propagation and 
constraint equations for this background, which describe 
the kinematics of the fluid flow and the dynamics of the 
gravitational field (see eqns. 87-100 in H3). These equa- 
tions are then split into their scalar, vector and tensor 
contributions giving three sets of linear propagation and 
constraint equations which govern the complete pertur- 
bation dynamics. 

From a dynamical systems point of view the past at- 
tractor Th for brane homogeneous cosmological models 
found in |9LllO| is a fixed point in the phase space of these 
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models. This phase space may be thought of as an in- 
variant submanifold within an higher dimensional phase 
space for more general inhomogeneous models. Our anal- 
ysis can be thought of as an exploration of the neighbour- 
hood of J~b out of the invariant submanifold explored 
in H03. 

We restrict our analysis to large - scales, at a time when 
physical scales of perturbations are much larger than the 
Hubble radius, A 3> H _1 (equivalent to neglecting Lapla- 
cian terms in the evolution equations). This may at first 
glance seem restrictive, but it is not the case for the non- 
inflationary perfect fluid models that are relevant to our 
discussion. Indeed, it is well known that any wavelength 
A < H^ 1 at a given time becomes much larger than H^ 1 
at early enough times. Because of this crucial property of 
perturbations for non - inflationary models our analysis is 
completely general, i.e. valid for any A as t — ► 0. 

In what follows we restrict our analysis to the case of 
vanishing background dark radiation (U = 0). 

DIMENSIONLESS VARIABLES AND 
HARMONICS 



quantities defined above: 

W: = icurlW *, S Qb = lcurlS ab , (5) 

8 a b = 77 curl 8 a b , TL a b — —H a b , (6) 
Id Id 

Ql = icurlC?:. (7) 

Finally, it is useful to use the dimensionless time deriva- 
tive d/dr = d/d [In a) (denoted by a prime) to analyse 
the past evolution of the perturbation dynamics. 

We use the harmonics defined in |22| to expand the 
above tensors X a , X a b in terms of scalar (S), vector (V) 
and tensor (T) harmonics Q. This yields a covariant 
and gauge invariant splitting of the evolution and con- 
straint equations for the above quantities into three sets 
of evolution and constraint equations for scalar, vector 
and tensor modes respectively: 

X = X S Q S (8) 
X a = k- 1 X s Q s a +X v Q v a (9) 
X ab = k- 2 X s Q s ab + k- 1 X v Q v ab + X T Q T ab . (10) 



We define dimensionless expansion normalised quanti- 
ties for the shear a a b, the vorticity uj a , the electric E a b 
and magnetic H a b parts of the Weyl tensor |2^| : 

r &ab TJ , k-"a c Eab n, Hab 

^ab — -Jj- , VVa = -j-jr , t a b = -Jp ' Hah = ~Jp ' 

where H is the Hubble parameter H = a/a and a is the 
scale factor. 

It turns out to be convenient to use the dimensionless 
variable 



w: 



aoj a 



(2) 



to characterise the vorticity of the fluid flow. Here cu a — 
VabcW bc is the usual vorticity vector. 

The appropriate dimensionless density and expansion 
gradients which describe the scalar and vector parts of 
density perturbations are given by (see j2f| for details of 
definitions) 



-B a p , Z* 



3n 
~H' 



(3) 



and for the brane - world contributions we define the fol- 
lowing dimensionless gradients describing inhomogeneity 
in the non - local quantities 



U* - 'flu n* - 



(4) 



where U a and Q a arc defined in equation (27) of 26] . 

When undertaking a complete analysis of the pertur- 
bation dynamics it turns out useful to define a set of aux- 
iliary variables corresponding to the curls of the standard 



PERTURBATION EQUATIONS AND THEIR 
SOLUTIONS 



We begin by giving the evolution equations for scalar 
perturbations (suppressing the label S): 



(11) 
(12) 
(13) 
(14) 
(15) 
(16) 



£' 


= (37- 


2)S- 


8, 




£' 


= (67- 


3)8- 


3 7 S, 




A' 


= (37- 


3)A- 


-7Z* , 




Z*' 


= (37- 


2)Z* 


-3(3 7 4 


•1)A- 


or' 


= (37- 


3)Q* 


-\U*- 


67A, 


u*' 


= (67- 


A)U* 







The scalar constraint equations are: 

2S = 2Z* - 3Q* , 
28 = 6A -3Q* + U* 



(17) 
(18) 



Vector perturbations on large - scale are described by the 
following evolution equations for the basic variables 

(19) 
(20) 
(21) 



(23) 
(24) 
(25) 
(26) 



£' 


= (37- 


2)S - 


8-2{j- 1)W* 


£' 


= (67- 


3)8- 


3 7 S + 67PF* , 


A' 


= (37- 


3)A- 




Z*' 


= (37- 


2)Z* - 


- (9 7 + 3)A - U 


Q*' 


= (37- 


3)0* 


- \U* -67A , 


U*' 


= (67- 


A)U* . 
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w*' 


= (37- 


A)W* 
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= (67- 


3)H- 


■ 12jW* + Q* , 
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TABLE I: Large-scale behaviour or geometric and kinematic quantities in the high energy limit. Here p — 37 — 3, q = 67 — 2, 

r — —3, s = 67 — 4, u = 67 — 5, v = 67 — 3 



Quantity Scalar contribution Vector Contribution Tensor Contribution 

A Al a« ,ig^a r .-^(g^y^a 5 A lfl ' ,4^I^V 

£T [7 *a s [7 *a s 

W W a u 

H Woa s Hoa s ,Hia2l 

w* 

Q* -W a s ,127W *a" 

E W a s Hoa s ,Hia2 q 

£ H a s H a s , -S-yHia^l 

HO 



and their curls 



w*' 


= (67 


-5)W* 
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(27) 


Q*' 


= (67 


-4)Q* 


- 367 2 W* , 


(28) 


n' 


= (97 


-4)W- 


6jW* 


(29) 


£' 


= (97 


-4)£- 


3 7 E, 


(30) 


E' 


= (67 


- 3)E - 


£ . 


(31) 



The constraints for the basic variables are 

3E = 4Z* - 6Q* + 6W* , (32) 

3£ = 12A - 6Q* + 217* , (33) 

H = 12-/W* - Q* , (34) 

E = H, (35) 

and the curl constraints are 

U = 6jW* , (36) 
E = 12 7 iy* - Q* , (37) 
£ = 12jW* - Q* . (38) 

Finally, the large-scale evolution of tensor perturba- 
tions are governed by propagation equations for the ten- 
sor contributions of the shear E a b, the electric £ a b and 



magnetic H a b parts of the Weyl tensor: 

E' = (3 7 - 2)E - £ , (39) 

£' = (67 - 3)£ - 3 7 E + H , (40) 

H' = (6-y-3)H-£, (41) 

together with propagation equations for their curls 

E' = (67 - 3)E - £ , (42) 

£' = (97 - A)£ - 3 7 E , (43) 

H' = (9 7 -4)7Y. (44) 



The only tensorial contribution to the constraints is 

E = H . (45) 

Solutions to these three sets of equations are presented 
in Table I. On close inspection of the exponents we con- 
clude that Th is stable in the past (as t — > 0) to generic in- 
homogeneous and anisotropic perturbations provided the 
matter is described by a non - inflationary perfect fluid 
with 7 - law equation of state parameter satisfying 7 > 1 
and we use the large-scale approximation (A 3> i? _1 ). 
In particular it can be easily seen that the expansion 
normalised shear vanishes as a — > 0, signalling isotropiza- 
tion. 

DISCUSSION OF RESULTS 

We have considered here only the case of vanishing 
background Weyl energy density, U = 0. This assump- 
tion considerably simplifies the analysis, but it is ex- 
pected that our results will remain true for U =/= 0. In- 
deed, when U ^ J-b still remains a past attractor of 
the isotropic models. In other words, our analysis is re- 
stricted to the invariant submanifold U = of the larger 
phase space with but this submanifold is asymp- 

totically stable against M / perturbations. A more 
complete analysis including this case will be the subject 
of a future investigation. 

In a related paper the dynamics of a class of spatially 
inhomogeneous G2 cosmological models in the brane- 
world scenario has been studied ■ Area expansion nor- 
malized scale - invariant dependent variables, the timclike 
area gauge and an effective logarithmic time t were em- 
ployed, and the initial singularity occurs for t — » —00. 
The resulting governing system of evolution equations 
of the spatially inhomogeneous G2 brane cosmological 
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models can then be written as a constrained system 
of autonomous first - order partial differential equations 
in two independent variables. The local dynamical be- 
haviour of this class of spatially inhomogcncous models 
close to the singularity was then studied numerically. It 
was found that the area expansion rate increases without 
bound (and hence the Hubble rate H — > oo as i ^ - oo, 
so that there always exists an initial singularity). For 
7 > 4/3, the numerics indicate isotropization towards Th 
as t — ► — oo for all initial conditions. In the case of radi- 
ation (7 = 4/3), the models were still found to isotropize 
as t — > —00, albeit slowly. From the numerical analy- 
sis we find that there is an initial isotropic singularity in 
all of these G2 spatially inhomogeneous brane cosmolo- 
gies for a range of parameter values which include the 
physically important cases of radiation and a scalar field 
source. The numerical results are supported by a qual- 
itative dynamical analysis and a detailed calculation of 
the past asymptotic decay rates |27j . 

Finally, we note that this result corrects an earlier pa- 
per [2^ in which two of the authors claimed the contrary 
to be tr ue ( i.e., brane -worlds do not isotropize). The 
claim in |22| was based on the existence of the decaying 
mode oc a~ 3 (growing in the past) in the solutions for 
A and Q* which when substituted into the shear Ijl7(l 
and divE constraints l|18|) gave the same mode in the ex- 
pansion normalised shear £ and electric part of the Weyl 
tensor £. This resulted from using an incorrect prop- 
agation equation for the non - local energy flux in |2j| 
(eqn. 32 in this paper) which gave different solutions for 
Q* [U- The error in |2(| results in a wrong decaying 
mode that is totally irrelevant for the future evolution 
of perturbations studied in that paper, but changes dra- 
matically their behaviour as t — » 0. It is easy to see that 
when the corrected solutions (in Table 1) are substituted 
into (|I7|) and i|18|) the a~ 3 mode actually cancels out in 
S and £. Note also that this mode does not appear in 
the solution for the magnetic part of the Weyl tensor Tt. 

Since S, £ and H give a better description of the geom- 
etry, their past evolution represents the true behaviour 
of anisotropics close to the initial singularity, and not A. 
In fact one can argue that the existence of the a~ 3 mode 
in A (defined in equation Q) results from it not being 
the most appropriate measure of inhomogeneity at high 
energies since the dominant background energy density 
is not p but p tot oc p 2 (see equation (11) in |29|). More 
precisely, if we define A^^ (High Energy) normalising 
with p tot instead of p then A HE oc A/p and the decaying 
mode in A becomes a mode oc a~ 3+37 . Since —3 + 37 > 
for 7 > 1, we remove the decaying mode in the density 
inhomogeneity. 
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